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ABSTRACT 

The purpose of this research was to show, both experiment- 
ally and analytically, that the stress concentration around 
holes in laminated fibrous composites is not adequately de- 
scribed by the theoretical solution for homogeneous orthotropic 
plates, and that it is a function of hole size. It was shown 
that the gross laminate properties, determined analytically from 
individual lamina properties, make it impossible to express the 
proper boundary conditions at a free edge. A new expression 
for the modulus of elasticity tangential to a free boundary 
was developed. Thus, the problem became similar to that of a 
circular hole strengthened by an elastic ring. Furthermore, 
the effective width of such a ring was found to depend upon 
hole size. Thus, the stress concentration factor was shown to 
be a function of hole size, which was not predicted by homogene- 
ous orthotropic theory. These results were very satisfactorily 
verified by experiment. 
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I. 



INTRODUCTION 



A. BACKGROUND 

Because of their high strength-to-weight ratios, laminated 
fibrous composites are expected to be utilized extensively for 
aerospace vehicles in the future. At the present time, how- 
ever, the realization of the maximum potential of composites 
is severely restricted by the limited techniques available for 
their analysis. An aerospace vehicle is necessarily a very 
complex structure which is difficult to analyze even when com- 
posed of homogeneous isotropic materials. The introduction of 
inhomogeneous anisotropic composites greatly magnifies the 
problem. Lekhnitskii presents a solution for stresses around 
elliptic holes in homogeneous orthotropic plates, but the ap- 
plicability to nonhomogeneous composites is questionable. 

B. OBJECTIVES 

The first objective of this research was to show experi- 
mentally that the homogeneous orthotropic solution for stress 
concentration around a circular hole is not accurate for 
laminated fibrous composites. The second objective was to 
determine whether hole size has an effect on the stress 
concentration in an infinite plate. The third objective was 
to analyze the effects of inhomogeneity on stress concentra- 
tion. Specifically, the boundary conditions on a free edge 
in an inhomogeneous laminate were analyzed in general and 
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in particular for the composite tested experimentally in this 
research . 

C . APPROACH 

As stated above, the experimental objectives required . 

the comparison of stress fields for different hole sizes with 

the homogeneous orthotropic solution. A (0°,-45°) s glass- 

epoxy composite was used. Four different hole sizes were 

tested by reflection photoelasticity. The analytical approach 

was to examine the relationships between constituent, lamina, 

2 

and laminate material properties given by Calcote . The pur- 
pose of this analysis was to determine if this attempt to re- 
present laminated fibrous composites as homogeneous materials 
was valid near free boundaries. 
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II. EXPERIMENTAL ANALYSIS 



A. SPECIMEN DESIGN 

1 . Material Properties 

The solution given by Lekhnitskii [3] for the stress 
concentration around a circular hole in an orthotropic plate 
is dependent upon the material properties. Therefore, the 
first task was to determine these properties experimentally. 

Due to the thickness of the photoelastic coating to be used 
for the stress concentration analysis, the entire glass-epoxy- 
photoelastic system was analyzed as a whole, rather than try- 
ing to apply correction factors. Thus, the material properties 
were also obtained by photoelasticity. 

Determination of the material properties of an ortho- 
tropic material requires tension specimens cut along the 
principle 1-t axes, and a third specimen cut at 45° to these 
axes [2]. Therefore, standard tensile specimens as described 
in the Advanced Composites Design Guide [1] were made up. 

2 . Stress Concentration 

Previous experimental work by this author seemed to 
indicate that stress concentration in composites is a func- 
tion of hole size. One purpose of the present research was 
to attempt to verify this hypothesis. Therefore, hole size 
was the only independent variable considered in designing 
specimens. All specimens were cut from the same sheet of 
material. The width of each specimen was twelve times the 
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hole radius. The dimensions of all specimens are given in 
Table II. A. 1. 



TABLE II. A. 1 

1200 SERIES SPECIMENS 

MATERIAL: GLASS/EPOXY/PSl-C 

G/E LAYUP: (0°, * 45° ) s 

THICKNESS: .15" 



PLATE # 


LENGTH 
L (in) 


WIDTH 
W (in) 


HOLE 

DIAMETER 

(2a) 


STRESS ANGLE 
(a) 


1200 


9 


1 


- 


0° 


1210 


9 


1 


- 


45° 


1220 


9 


1 


- 


90° 


1203 


18 


3 


1/2 


0° 


1205 


12 


1-1/2 


1/4 


0° 


1206 


7-1/2 


3/4 


1/8 


0° 



B. SPECIMEN PREPARATION 

The G/E composite was obtained from 3-M Company in a two 
by three foot sheet. The specimens were cut from this sheet 
using a diamond edge saw. All cutting operations were water 
cooled. The ends of each specimen were built up by attaching 
pieces of (0°, 90°) G/E as described in the Composites Design 
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Guide. The most difficult operation was that of drilling 



holes in the G/E material. The difficulty lay in preventing 
fiber pullout as the drill bit came out the back of the 
material. This problem was overcome by sandwiching the 
composite between aluminum when drilling. 

The PS-1C photoelastic coating was cut one-eighth inch 
oversize and the holes were drilled. After cementing the 
coating to both sides of the G/E composite, the holes were 
reamed and the edges were trimmed by machining at high speed. 

C. TESTING PROCEDURE 

The photoelastic data was taken using the 030 series re- 
flection polariscope by Photolast ic, Inc. The specimens 
were mounted in a Reihle 10,000 lb. testing machine. In 
order to reduce the scatter due to viscoelastic effects, 
specimens 1200, 1210, and 1220 were cycled to 300 lb. three 
times before testing began. The other specimens were held 
at their test load for five minutes before taking data. The 
tare values were taken after this initial stretching. 

D. MATERIAL PROPERTIES RESULTS 

The original and reduced data are given in tabular form 
in the Appendix. Note that the fringe value, F, has been 
divided by the scale factor of 47 of the linear compensator. 
Thus, the "fringe orders" given (i.e. NN, NO, TNN, and TNO) 
are actually compensator readings. Using these data, the 
stress-strain curves for each specimen and the Poisson's 
Ratio vs. strain curves for specimens 1200 and 1220 were 
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point plotted (Figures II. D. 1-5). There is considerable data 
scatter, particularly in Poisson's Ratio. Although a good deal 
of effort was directed towards eliminating this scatter', visco- 
elastic effects and slight variations in load during testing 
made the task extremely difficult. It appears that the visco- 
elastic properties of the material affect the transverse 
strain to an even greater degree than the longitudinal strain. 
In view of these difficulties, due consideration was given to 
the relationship. 



E 1 v 21 = E 2 v 12 (II.D.l) 

when the curves were drawn in these figures. Using the equa- 
tion 




(II. D. 2) 



the following material properties were derived from the 
curves : 



E 1 “ !* 24 MSI v 12 = >5 

E 2 = °- 78 MSI v 21 = . 31 

G 12 = .51 MSI 

E. STRESS CONCENTRATION RESULTS 

A point-by-point photoelastic analysis was done for plates 
1203, 1205, and 1206. The points chosen were along the X and 
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FIGURE II.D.l. 
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FIGURE II. D. 2 - Modulus at 45 Degrees 
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FIGURE II. D. 3 - TRANSVERSE MODULUS 
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FIGURE II. D. 4 - Longitudinal Poisson's Ratio 
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FIGURE II. D. 5 - Transverse Poisson's Ratio 
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Y axes as shown in Figure II.E.l. The data are tabulated in 
the Appendix as explained in the previous paragraph. The 
longitudinal strain variation along the X-axis from the hole 
to the edge of the plate is plotted in Figures II. E. 2-5. ■ 

The theoretical strain concentration at the hole is given 
by Lekhnitskii as 

K = p (1+n) (II.E.l) 



where 



n 




v 12> 




(II. E. 2) 



These equations give a theoretical strain concentration factor 
of 3.15 for this material. It is evident from Figure II. D. 5 
that the strain concentration is a function of hole size for 
this laminated fibrous composite, even though the hole-size- 
to-plate-width ratio is the same. Also, the measured value 
seems to approach the theoretical value as hole size increases. 
Previous tests by this author on plates of constant width gave 
the same results. 

Figure II. D. 5 also shows that the higher the strain con- 
centration at the hole, the faster it drops off with the 
normalized distance from the hole. In fact, the former is 
actually lower at some points! This phenomenon is necessary 
so that the average value of the stress concentration over 
the entire X-axis (including the hole) will be equal to one. 
Note also the drop in strain concentration at the outside edge 
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of each plate. This effect is due to a change in gross 
laminate properties at a free edge, which is discussed in de- 
tail in Section III. 
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III. THEORETICAL ANALYSIS 



The solution given by Lekhnitskii for the stresses around 
elliptic holes in orthotropic plates assumes a homogeneous 
material. Because the stress distribution is a function of 
the material properties, it is necessary to formulate gross 
laminate properties which approximate the elastic properties 
of the composite by those of a homogeneous system. The obvi- 
ous question here is, "Are the basic assumptions necessary 
for this formulation valid throughout the plate, particularly 
near a free boundary?" 

A. BASIC FORMULATION 

The basic formulation will first be derived in detail be- 
cause it is the basis for the discussion which follows. The 
stress-strain relationship (referred to natural axes) for an 
orthotropic lamina is 

{a} = [C] {e } (III.A.l) 

where 




and 
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[C] = 



C 11 C 12 0 



C 12 C 22 0 






'44 



Referred to arbitrary axes, this relationship becomes [2] 



r°r 




_ i 


] 


\°2 


- 


- ic] ; 


e 2 


L T 12 

a 




1 





(III. A. 2) 



where 



[C] = [T] —1 [C] [T] 



and [T] is the matrix of direction cosines; 





2 

cos a 


. 2 
sin a 


2sina 


cosa 


[T] = 


. 2 
sin a 


2 

cos a 


-2sina 


cosa 




-since cosa 


sina cosa 


2 

cos a 


. 2 
-sin 



(III .A. 3) 



Equilibrium requires that (assuming unit width) 

{P } C = £{P} (k) = £ [a] (k) h (k) 
where h^ k ^ is the thickness of the lamina. Thus, 



„ ,(k). (k) 

{a} c = liaj h 



Eh 



(k) 



(III. A. 4) 



Substituting equation III. A. 2 into this expression yields 

[C] C {e} C = Z [ ^ ] (k) (III. A. 5) 

Eh 
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At this point, the assumption is made that the strain field 
in each lamina is identical to that of the composite as a 
whole, that is 



{e} = {e} c = {e} (k) (III. A. 6) 

equals a constant over the summation. Thus, equation III. A. 5 
reduces to 



„ v rri 00 u (k) 

[C ] c = LA£1 h 



Eh 



(k) 



(III. A. 7) 



B. FREE EDGE BOUNDARY CONDITIONS 

The boundary condition for a free edge in a homogeneous 
plate loaded in-plane is 

a = 0 (III.B.l) 

n 



By applying this boundary condition to a composite represented 
by the homogeneous properties given by III. A. 7, one is actually 
only requiring that the average stress across the thickness be 
zero, rather than the stress at each point across the thick- 
ness of the plate. Applying equation III. A. 4 to III.B.l, 



c Ea (k) h (k > 

° = — TFT— = 0 

n Eh k 



(III.B.2) 



(k ) ' 

Thus, a n in general need not be zero, whereas the physics 

of the problem require that it must be zero. Therefore, the 

material properties given by equation III. A. 7 are not valid 

near a free boundary in a laminated composite. 
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Equation III.B.2 results in the free-edge displacement 
shown in Figure III.B.2. The free edge is still a plane nor- 
mal to the surface. Thus, all strains and displacements are 

(k) 

constant through the thickness. The stress a n in a lamina 
is derived as follows. The tangential stress, a^, is known 
and the normal stress, a n , is zero. Thus, at the free edge 



(k) c c ,„c 

e = £ = a / E 

s s s s 



(III.B.3) 



and 



(k) 



n 



= £ 



n 



c c 
-v £ 
sn s 



But 



(k) _ 



n 



= ( 



n 



1 - v v 
ns sn 



-) (k) (£ c +v (k) £ C ) 
n sn s 



or 



00 _ 



n 



E 



= ( 



n 



1 -v v 
ns sn 



. . (k) c . c 

-) (v -v ) £ 

sn sn s 



(III.B.4) 



Thus, the error in the boundary condition is directly propor- 
tional to the difference in Poisson's ratios between the 
composite and the lamina in question. The error is zero for 
homogeneous materials. 

The free edge actually looks like Figure III.B. 3 . The 
displacement must be continuous between lamina. Thus, the 
displacement will vary through the thickness of a given 
lamina . 
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m.B.l - Coordinate System 




Fig. III.b. 2 - Displacement for a C = o 

n 




Fig. III.b. 3 - Displacement for 0 

n 
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C. MATERIAL PROPERTIES AT A FREE BOUNDARY 

It has been shown that equation III.B.2 is not the correct 
boundary condition for a free boundary in a laminated compos- 
ite. The correct condition is 

o (k) = 0 for all k (III.C.l) 



In order to derive the composite material properties, an addi- 
tional assumption is required: 




(III.C.2) 



is constant through the thickness. Note that 



e 



(k) _ 
n 



'n 



Hooke's Law for this case is 



a 



c 

s 



E c e 
s s 



and 



a 



(k) _ 

s 




(III.C. 3) 



(III.C.4) 



(III.C. 5) 



Recalling equation III. A. 4 



So (k, h (k > 
c _ _s 

° s “ 



Eh 



or 



E C e 
s s 



EE (k) e h (k) 
s s 



Eh 



(k) 
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Thus 



E 



c 

s 



EE (k) h (k) 

s 



Zh 



(k) 



(wc.&) 



D. NUMERICAL RESULTS 

The next step in the analysis was to find numerical values 

c c 

for E s , Young’s modulus at a free boundary and Eg, Young’s 

modulus in an arbitrary direction far from any boundaries. 

Q 

E g may be found by use of equation III.C.6, which requires that 
E g be known for each lamina. A single lamina composite 
could have been made up and tested, but the data scatter in- 
herent in the properties of glass-epoxy composites cut from 
different sheets would have rendered the results virtually 
useless for this analysis. Fortunately, this procedure was 
not necessary. It was possible to analytically reduce the 
gross laminate properties to individual lamina properties by 
the procedure which follows. The experimentally determined 
properties of the G/E/PS1-C composite are repeated here for 
convenience . 

E^ =1.24 MSI vj 2 = .5 

E 2 = .78 MSI V 21 = .31 

G^ = .51 MSI 

The stiffness matrix may be computed from 
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(III.D.l) 



'11 



1_V 12 V 21 



'22 



E 2 

1 " V 12 V 21 



E 2 V 12 



'12 



1-v 12 V 21 



C 44 G 12 



which give 



[C] = 



1.47 

.462 

0 



.462 

.923 



0 

0 



.51 



(III.D.2) 



The properties of the PS1-C photoelastic coating are 
E P = .462 



v p = .36 



and the stiffness matrix is 
.528 .190 

.190 .528 

0 0 



(C] P = 



0 

0 



.169 



(III.D.3) 



The thickness of the G/E and PS1-C portions of the composite 
will also be required. 

h C = .15 

h C-p = .0625 



h P = .0875 



The thickness of the PS1-C includes the cement, which has ap- 
proximately the same properties. Recalling equation III. A. 7, 
it is possible to determine the properties of the G/E system 
without the photoelastic coating. 

h C [C] C = h P [C] P +h C-p [C] C-p (III.D.4) 



38 



yields 



[C] 



c-p _ 



2.79 
. 843 

0 



.843 

1.48 

0 



0 

0 



987 



The inverse of equations III.D.l are [2] 

„2 



E 1 C 11 C 



12 



E 2 C 22 “ C 



22 

2 

12 



11 



(III.D.5) 



v 



12 



12 

'22 



C 11 



(III.D.6) 



G 12 C 44 



which gives the properties of the G/E system. 



E C " P = 
E 1 


2.31 


MSI 


v c-p = 
12 


.57 


M 

ro 

0 

1 

"0 

II 


1.23 


MSI 


c-p 

V 21 


. 30 



G ^ P = .987 MSI 



Recalling equation III. A. 7 and noting that the thickness of 
each lamina in the G/E system is the same: 



[C] 



c = Z [C] (k) 
n 



(III.D.7) 



where n is the number of lamina. 
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Recalling that 



[C] = [T] 1 [C] [T] (III.D.8) 
it is possible to solve for [C] because all lamina are the 
same except for orientation! Equation III.D.8 may also be 
written [2] 



- 4 2 2 4 

C^=C^cos a+2 (C^ 2 + 2 C^^) sin acos a+C 22 s; *" n a 

- 2 2 4 4 

C^ 2 = (C^^+C22“4C^^ ) sin acos a+C^tsi- 11 a+cos a) 



- 4 2 2 2 

C' 22 =C ll S i' n a+2 (C^2 + 2C^^ ) sin acos a+C 22 COS a (III.D.9) 



_ 2 2 4 4 

(C^^+C 22 -2C jL 2 -2C 44 ^ s; *- n acos (sin a+cos a) 



Noting that the C^'s do not vary across the summation in 
equation III.D.7, one may proceed as follows. The orientation 
of the various lamina are 



c' 1 ' = a (6) = 0* 





(2) 

a 


- cx< 5 > = 


+ 45° 








(3) 

a 


- a W = 


-45° 






and 


equation III. 


D.7 reduces 


to 




M 


C-p 


“3 


2 


1 


4 




C 12 




1 


4 


1 


-4 


-< 


C 22 


► = 1/6 


1 


2 


"3 


4 




c, , 

v 44 

^ -J 




1 


-2 


1 


2 J 




-\ 

11 



12 



V- 



22 



44 



(III.D.8) 
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The inverse is 




(III.D.9) 



which gives the stiffness matrix for a single lamina as 



[C] = 



4.78 

.161 

0 



.161 

.852 

0 



0 

0 

.305 



(III. D. 10) 



Applying equation III.D.6 gives the material properties of 
the lamina . 

=4.75 MSI \> lt = .19 



e t = . 847 MSI v t£ = .03 

G £t = .305 MSI 

(k) 

In order to apply equation III.C.6, E must be known for 

s 

each lamina. This quantity may be computed from 



e Z 20 , C Z . 20 , , r c Z 0 .. -j . 2 O 0 

— Tk) = cos + 77 Sln + 1/4 >57 • 2v it lsln 2 

C t Aj t 



or 



1 . 211cos 20 + 1 . 18sin 20 + .8sin 2 (2 0 ) (III. D. 11) 



' £ 



(k) 
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01*5 




GGC 0C2 QCi 0C5 039 



ANGLE, 0 



K-SCRLE=2. 00E+01 UNITS INCH. 
V-5CRLE--5. 00E-01 UNITS INCH. 



FIGURE III.D.l - Young's Modulus Curves 

SRBR0072 
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( k ) G 

E is computed for each lamina and then E is computed for 

S S 

equation III.C.6. 

Q 

For comparison, E may be computed from 



- = cos 2 9 + 1.59 sin 2 9 + .358 sin 2 (20) 



or 



.806 cos 2 9 + 1.28 sin 2 9 + .289 sin 2 (29) 



(III. D. 12) 



c c 

Eg and E g are plotted as a function of 0 in Figure III.D.l. 



E. EFFECT ON STRESS CONCENTRATION 

As seen in Figure III.D.l, Young's modulus in the y direc- 
tion at the boundary is 1.17 MSI compared to 1.24 MSI far from 
all boundaries. The variation along the X-axis is something 
like Figure III.E.l. Thus, the material at the boundaries is 
less stiff than the material near the center of the plate. 

Now, suppose that the theoretical stress concentration varia- 
tion for constant material properties is given by curve A in 
Figure III.E.2. The variation in material properties described 
by Figure III.E.l would tend to reduce the stress concentra- 
tion near the boundaries, because the material is less stiff 
in these areas (Curve B) . This effect is shown clearly at 
the boundaries of the plates in Figures II. D. 2-5. Thus, the 
lower stiffness near a free boundary in a laminated composite 
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requires that the stress concentration factor be smaller than 
that given by homogeneous orthotropic theory. 

The effect described above can also be shown to vary with 
the hole size. The modulus of elasticity of the laminated 
composite is lowest at the free boundary and increases to a 
steady value at a finite distance from the boundary as shown 
in Figure III.E.l. The actual function will depend on compon- 
ent materials and layup of composite, but will depend only on 
the X distance from the free edge for a particular composite: 

E 1 = E ± (X) (III.E.l) 

The stress concentration along the X-axis, however, is a 
function of the number of hole radii away from the boundary: 

K = K(X/r) (III.E.2) 

Thus, in terms of the normalized distance, X/r, a plate with 
a small hole is affected more by the difference in stiffness 
than a plate with a large hole. Therefore, the smaller the 
hole, the smaller the stress concentration factor, even in an 
infinite plate or plates with a constant hole-size-to-plate- 
width ratio. 

If Young's modulus in Figure III.E.l is assumed to vary 
like Curve B, rather than Curve A, one has the problem of 
stress concentration around holes which are strengthened by 
elastic rings. This problem has been solved by Savin [4] who 
describes similar results. The theoretical stress concentra- 
tion for infinite copper plates with holes strengthened by 
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FIGURE III. E. 1 - E near a free boundary 




FIGURE III.E.2 - Inhomogeneous effects on k 
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steel rings of various widths is presented by Savin and this 
set of curves is very similar to Figure II. E. 5, the experi- 
mental results for plates with different size holes. Consider 
a series of such rings, one inside the other. A large number 
of thin rings with appropriate material properties could be 
used to produce the variation in stiffness described by Curve 
A of Figure III.E.l. It is suggested that this approach be 
further investigated in order to obtain numerical results. 
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IV. CONCLUSION 



The first objective of this research was to show experi- 
mentally that the homogeneous orthotropic solution for stress 
concentration around a circular hole is not accurate for 
laminated fibrous composites. Referring to Figure III.E.5, 
the following stress concentration factors are obtained: 



TABLE IV. 1 



Theory 3.15 
Plate 1203 (r = .25") 3.07 
Plate 1205 (r = .125") 2.45 
Plate 1206 (r = .0625") 2.01 



Obviously, the theory is not accurate. The second objective 
was to determine whether hole size has an effect. on the stress 
concentration in an infinite plate. Referring to Table IV. 1 
again, it is obvious that the stress concentration factor de- 
creases with hole size. 

The third objective was to analyze the effects of in- 
homogeneity on stress concentration. It was found that to en- 
force the boundary condition at a free edge required the 
modification of the gross laminate properties near the edge. 
The lower stiffness obtained analytically required that the 
stress concentration factor be lower than that predicted by 
homogeneous orthotropic theory. Furthermore, it was shown 
that the stress concentration factor varies with hole size 
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because the modified material properties affect a larger 
normalized area for smaller holes. 
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APPENDIX A: EXPERIMENTAL DATA 



A- I MATERIAL PROPERTIES TEST DATA 
MATERIAL PROPERTIES TEST 9 APR 



PLATE # 1200 



# OF 


RUNS: 4 


AREA 


0.150 


F 


= 40 . 


248 






P 


S I GMA 


NN 


NO 


OELNN 


DELNO 


EPSY 


EPSX 




0 


0 


1 . 1 


-2.7 












100 


666 


22.9 


2 0.2 


21.8 


22.9 


505 


-3 72 


* 


200 


1333 


43.0 


46.6 


41.9 


49. 3 


1239 


-396 


• 


300 


2000 


62.5 


66. 6 


61.4 


69.3 


1712 


-758 


• 


40 0 


2 666 


78.7 


84.4 


77.6 


87.1 


2135 


-988 


• 


500 


3333 


94.5 


106.1 


93. 4 


108.8 


2309 


-949 


• 


0 


0 


0.4 


-2.7 












100 


666 


22.6 


19.7 


22.2 


22.4 


458 


-434 


• 


200 


1333 


41 . 5 


43.2 


41 . 1 


45 . 9 


1116 


-52 7 


• 


300 


2000 


61.4 


64.9 


61.0 


67.6 


1626 


-829 


• 


400 


2666 


79. 8 


85.6 


79.4 


89.3 


2195 


-1000 


• 


500 


3333 


96.7 


107.5 


96.3 


110.2 


2777 


-1098 


c 


0 


0 


-1 .4 


-3.0 












100 


666 


20.9 


21 .0 


22.3 


24.0 


55 1 


-346 


• 


200 


1333 


40.2 


42.0 


41.6 


45.0 


1042 


-6 31 


• 


300 


2000 


61. 1 


65.4 


62.5 


68.4 


1613 


-901 


• 


400 


2666 


81.2 


90.4 


82.6 


93. 4 


2314 


- 1 0 1 0 


%> 


500 


3333 


98. 1 


106.4 


99.5 


109.4 


2600 


-1404 


• 


0 


0 


2.0 


1 . 7 












100 


666 


19.6 


10.2 


1 7.6 


8.5 


-195 


-903 


•A. 


200 


1 333 


42.5 


45.0 


40.5 


43.3 


984 


-645 


• 


300 


2000 


60.3 


66.5 


53.8 


64.8 


1545 


-321 


• 


400 


2666 


80. 1 


87.6 


78.1 


85.9 


2042 


-1100 


• 


500 


3333 


99. 1 


109.3 


97.1 


107.6 


2587 


-1320 


• 



75 



NU 



737 

307 

443 

463 

333 

948 
481 
5 10 
4 55 
395 

628 

606 

559 

436 

540 

O,. J. J- 

'f 'n v 

655 

531 

539 

51D 
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MATERIAL PROPERTIES TEST 



6 MAY 75 



PLATE ft 1210 



# OF 


RUNS: 4 


AREA 


0.150 




= 40. 


248 






P 


SIGMA 


NN 


NO 


DELNN 


DELNO 


EPSY 


EPSX 


NU 


0 


0 


-0.3 


-0 . 8 












100 


666 


17.1 


19.4 


17.4 


20.2 


519 


-181 


. 349 


200 


1333 


33.0 


38.3 


33.3 


39. 1 


1020 


-319 


.313 


300 


2000 


49.0 


57.3 


49.3 


58.1 


1523 


-460 


.30? 


400 


2666 


66.0 


77.2 


66.3 


78.0 


2040 


-627 


.307 


500 


33^3 


81.7 


94.3 


82.0 


95. 1 


244 1 


-859 


.352 


0 


0 


1 . 4 


-3 . 5 












100 


666 


18.2 


20.4 


16.8 


23.9 


766 


90 


OU -JU- a- 
I 1 'T -r 


200 


1333 


34. 1 


40.5 


32.7 


44.0 


1340 


24 


«£■ tX, JU 

t 'i' -r r 


300 


2000 


51 .0 


60.4 


49.6 


63.9 


186 1 


-134 


.072 


400 


2666 


66. 0 


81.0 


64.6 


34.5 


250 1 


-98 


.039 


500 


3333 


81 .6 


98.2 


80.2 


101.7 


291 1 


-315 


.108 


0 


0 


-1.3 


-1 .0 












100 


666 


20. 8 


23.5 


22. 1 


24.5 


539 


- ?c;cj 


.508 


200 


1333 


33.4 


39.4 


3 9.7 


40.4 


841 


-756 


. 899 


30 0 


2000 


53.4 


5 6.5 


54. 7 


57.5 


1269 


-931 


. 734 


400 


2666 


71.0 


77.0 


72.3 


73.0 


1799 


-1 1 10 


.617 


500 


3333 


83 . 8 


97.0 


85.1 


98.0 


2491 


-933 


.375 


0 


0 


-1.3 


-1 .0 












100 


666 


24.5 


19.5 


25.8 


20.5 


199 


-839 


^ nO 
't' 't" +C* V 


200 


1333 


37.7 


41 .0 


39.0 


42.0 


955 


-603 


.625 


30 0 


2000 


55.3 


58.7 


5 6.6 


59.7 


1326 


-951 


.717 


40 0 


2666 


71.0 


80.1 


72.3 


81. i 


1986 


-92 3 


.465 


500 


3333 


87.0 


95 . 1 


88. 3 


96. 1 


2247 


-1306 


.5 81 
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MATERIAL PROPERTIES TEST 
PLATE # 1220 



6 MAY 75 



# OF 


RUNS: 


4 AREA 


0.150 


F 


• 

o 

Nj“ 

II 


248 






P 


SIGMA 


NN 


NO 


DELNN 


OELNO 


EPSY 


EPSX 


NU 


0 


0 


1 .0 


-2.0 












100 


666 


29. 8 


28.2 


28.8 


30.2 


564 


-495 


. 745 


200 


1 333 


55.9 


65.2 


54.9 


67.2 


1847 


-362 


. 1 95 


300 


2000 


82 . 7 


94.6 


81 .7 


95.6 


2543 


-744 


.2 9° 


AO 0 


2666 


1 09.0 


127.2 


108.0 


129.2 


3453 


-893 


.259 


0 


0 


2.6 


1 .0 












100 


666 


31.0 


29.8 


23.4 


2 3.8 


595 


-547 


.919 


200 


1 333 


55.4 


66.4 


52.8 


6 5.4 


1823 


-301 


. 1 65 


300 


2000 


33.8 


96.8 


81.2 


95.8 


2515 


-752 


.2 99 


400 


2666 


1 10.0 


130.2 


107.4 


129.2 


3477 


-845 


.243 


0 


0 


3 . 4 


-1.2 












100 


666 


32 . 5 


31.7 


29.1 


32.9 


815 


-356 


.437 


200 


1333 


59. 0 


6 4.6 


55.6 


65.8 


1734 


-503 


. 290 


300 


2000 


86.7 


98.3 


83.3 


99.5 


2654 


-693 


.2 63 


400 


2666 


112.8 


128.6 


109.4 


129.8 


3433 


-969 


.282 


0 


0 


2.6 


-1 .6 












100 


666 


33.5 


28 . 2 


30.9 


29.8 


555 


-688 


A vU •'* 

*Y* -i' 'r 


200 


1 333 


56. 4 


64.5 


53.8 


55 . 1 


1825 


-340 


. 1 86 


300 


2000 


86.6 


95.0 


3 4.0 


97.6 


251 1 


-369 


. 346 


400 


2666 


112.7 


129.0 


110.1 


130.6 


3453 


- 978 


.2 83 
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A-2 STRESS CONCENTRAT I CN TEST DATA 



STRESS CONCENTRATION DATA 

PLATE n 12D3 X-AXIS 

# OF RUNS: 2 LOAD 800.0 f-= 36.596 



29 MAY 75 



X/P 


TNN 


T\'0 


NN 


NO 


EPSY 


EPSX 


K 


A 


FF 


- 2.0 


-6.5 


43.5 


43.8 


1096 


-569 


1 .00 


0.52 


1.0 


-8.0 


-6.5 


120.4 


jfc # >;c # s}c 3J; 


3455 


-v * 

*1 % -1' 'Y' "V '(• 'i' 


3.15 


•J+ n*- J/ 

r 


1.3 


-2.0 


-6.5 


75.0 


75.4 


1677 


-1 139 


1 . 53 


0.68 


1 . 5 


-2.0 


-6.5 


55.5 


60.5 


1573 


-530 


1.44 


0. 34 


1.8 


-2.0 


-6.5 


45.6 


49.6 


1337 


-404 


1 . 22 


0.30 


2.0 


-2.0 


-6.5 


45.6 


46.4 


1161 


-580 


1 .06 


0.50 


2.5 


- 2.0 


-6.5 


45.6 


46 . 4 


1161 


-580 


1.06 


0.50 


3.0 


-2.0 


-6.5 


A 5 . 6 


46.4 


1161 


-580 


1.06 


0.50 


4.0 


-2.0 


-6.5 


45.6 


46.4 


1161 


-580 


1.06 


0.50 


5.0 


-2.0 


-6.5 


45.6 


46.4 


1161 


^ j580 


1.06 


0.50 


6.0 


-8.0 


-6.5 


30.7 


.V -'t J. 

•V -v 'i* 'A T* *>' 


1041 


^ 


0.9 5 


a- j, 

V V V ' 


FF 


-2.0 


-6.5 


43.5 


4 3.8 


1 096 


-569 


1.00 


0.52 


1 .0 


-5.8 


-6.5 


113.2 


U« ^ V/ vW ^4 

-I- "l- '1' '1- 'I* 


3202 


'l* *! X V '1* '» x 


2.92 


*•< »V X J. <» «<. 
•V" 'S' 'i' -1' 'l' 


1.3 


-2.0 


-6.5 


75.7 


87.0 


2 2 52 


-627 


2.0 5 


0.28 


1.5 


-2.0 


-6.5 


56.8 


73. 1 


2217 


65 


2.02 


-0. 02 


1.8 


-2.0 


- 6.5 


47.3 


57.7 


1 720 


-34 


1.57 


0.05 


2.0 


-2.0 


-6.5 


47.3 


50.0 


1 297 


-506 


1.18 


0.39 


2.5 


-2.0 


-6.5 


47.3 


50.0 


1297 


-506 


1.18 


0.39 


3.0 


-2.0 


-6.5 


47.3 


50.0 


1 297 


-506 


1.13 


0.3 9 


4.0 


-2.0 


-6.5 


47. 3 


50.0 


1 297 


-506 


1.18 


0. 39 


5.0 


-2.0 


-6.5 


47.3 


50.0 


l 297 


-506 


1.18 


0.39 


6. 0 


- 8.0 


-6.5 


3 5.4 


sjc ff. if y’c 


1167 


,1. X X 

v 'i' -V- -v n' 


1 .05 


x j. x x. x x 

V *1' 'X V *>' 
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STRESS CONCENTRATION DATA 



28 MAY 75 



PLATE # 1205 X-AXIS 



# OF 


RUNS: 


2 


LOAD 4 


00.0 


o 

st 

II 

u. 


.248 






X/R 


TNN 


TNO 


NN 


NG 


EPSY 


EPSX 


K 




FF 


-0.5 


-5.0 


51.7 


53.4 


1 424 


-676 


1.00 


0.47 


1.0 


-2.5 


-5.0 


117. 5 


u. wt. .U J. X 
T“ <V —l* 


3551 


^ a j, ^ j, 

'r nr ' r t r f 


2.49 


J. X J/ ^ 

•v or 5-* x -r 


1.5 


-0. 5 


-5.0 


63.2 


82.1 


2694 


130 


1.89 


-0.05 


2.0 


-0.5 


-5.0 


49.9 


51 . 7 


1 394 


-633 


0.98 


0.45 


2.5 


-0.5 


-5.0 


49.9 


51.7 


1 394 


-633 


0.98 


0.45 


3.0 


-0. 5 


-5.0 


49.9 


51.7 


1394 


-633 


0.98 


0.45 


4. 0 


-0.5 


-5.0 


49.9 


51.7 


13 94 


-633 


0.98 


0.45 


5.0 


-0.5 


-5.0 


49.9 


51.7 


1394 


-633 


0 . °8 


0. 45 


6. 0 


-6. 5 


-5.0 


36.1 


X X X X X 
^ "r v 'f* 


1260 


X X X X %•/ 

v -r* 


0.88 


X X X X X 

v 


FF 


-0. 5 


-5.0 


51.7 


53.4 


1424 


-676 


1.00 


0.47 


1.0 


-5.2 


-5.0 


111.0 


J* jly jC Jv 


3438 




2.41 


»•» <x> 

nr 'r -r 'r x -r 


1.5 


-5.2 


-5.0 


50.3 


76^0 


2636 


382 


1.85 


-0. 14 


2.0 


-5.2 


-5.0 


49. 1 


51.0 


1195 


-990 


0.84 


0.83 


2.5 


-5.2 


-5.0 


49. 1 


51.0 


1 195 


-990 


0.34 


0.8 3 


3.0 


-5.2 


-5.0 


49.1 


51.0 


1 195 


-990 


0 . 84 


0.83 


4.0 


-5.2 


-5.0 


49. 1 


51.0 


1 195 


-990 


0.84 


0.83 


5.0 


-5.2 


-5.0 


49.1 


51.0 


1195 


^-990 


0.84 


0.83 


6. 0 


-2.6 


-5.0 


40. 2 


»V jV ^ 

'l‘ -V -V v- *ir -V* 


1266 


^ ^ ^ ^ V 


0.89 


X X J# v*/ X X 

or* -*r -V* 'i' -*r 
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STRESS CONCENTRATION DATA 



29 MAY 75 



PLATE # 1206 X-AXIS 

4 OF RUNS: 2 LOAD 200.0 



X/R 


TNN 


TNC 


NN 


NO 


FF 


- 1.0 


-5.5 


55.6 


61. 


I .0 


-8.0 


-5.5 


109.0 


^ V 


2.0 


-1.0 


-5.5 


46 . 6 


67. 


3.0 


-1.0 


-5.5 


46.6 


57. 


4.0 


-1.0 


-5.5 


46.6 


57. 


5.0 


-1.0 


-5.5 


46.6 


57. 


6.0 


-5.8 


-5.5 


40. 1 


•J* V* 

T* *V> *V -V •¥“ 


FF 


-1.0 


-5.5 


55.6 


61. 


1.0 


-8.0 


-5.5 


112.0 


*4-* X X 

X* 


2.0 


-1.0 


-5.5 


52.7 


64. 


3.0 


-1.0 


-5.5 


52.7 


53 . 


4.0 


-1.0 


-5.5 


52.7 


53. 


5.0 


-1.0 


-5.5 


52.7 


53. 


6.0 


- 5. 8 


-5.5 


41.6 


4^ 4# 

V -v 'r 



F = 36.596 



EPSY 


EPSX 


K 


A 


15 90 


-481 


1.00 


0.30 


3 148 


r ^ ^ V r r 


1 . 98 


X X X 4/ X 

n 1 'r v “r X* V 


2265 


523 


1.42 


-0.23 


1710 


-31 


1 .08 


0. 02 


1 710 


-31 


1.08 


0. 02 


1710 


-31 


1.08 


0.02 


1235 


4# X X X 
^ 


0.78 


£ * # # * 


1 590 


-431 


1 .00 


0.30 


3229 


X X X X X J/ 
X* T ' ' 'J' 


2.03 


OL. JL. vU vC- 

'O Y 'I' Y Y ** 


1855 


-109 


1.17 


0.0 6 


1273 


-6 91 


0.80 


0. 54 


1273 


-691 


0 .80 


0.5 4 


12 73 


-691 


0.30 


0. 54 


1275 


v V V v ^ 


0.80 


X X X X *v x 
; r t nr* 



2 

5 

4 

4 

4 

2 

T- 

1 

5 

5 

5 

.v 

'r* 
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STRESS CONCENTRATION DATA 



29 MAY 75 



PLATE # 1203 Y-AXIS 

# OF RUNS: 2 LOAD 300.0 



Y/P 


TNN 


TNG 


NN 


NO 


FF 


-2.0 


-6.5 


43.5 


43. 


1.0 


15. 5 


-6.5 


68.7 


'r' 'r '7 -r v 


1.3 


-2.0 


-6.5 


20.4 


7. 


1 . 5 


-2.0 


-6.5 


14. 3 


4. 


1.8 


-2.0 


-6.5 


10.2 


9. 


2.0 


-2.0 


-6.5 


16.7 


14. 


2.5 


-2.0 


-6.5 


26.4 


22. 


3.0 


-2.0 


-6.5 


30.1 


30. 


4.0 


-2.0 


-6.5 


38.0 


33. 


5.0 


-2.0 


-6.5 


40.0 


39. 


6.0 


-2.0 


-6.5 


40.0 


39 . 


FF 


-2.0 


-6.5 


43.5 


43 . 


1.0 


14.0 


-6.5 


70.5 


'V ■**r OU 
'T T T T* 'i' 


1.3 


-2.0 


-6.5 


24.3 


13. 


1.5 


-2.0 


-6. 5 


14.4 


2. 


1 . 8 


-2.0 


-6.5 


14. 6 


6. 


2.0 


-2.0 


-6.5 


22.1 


14. 


2. 5 


-2.0 


-6.5 


26.9 


22. 


3.0 


-2.0 


-6.5 


33.2 


32. 


4.0 


- 2.0 


-6.5 


38.2 


38. 


5.0 


-2.0 


-6.5 


42.6 


41 . 


6. 0 


-2.0 


-6.5 


44.5 


43. 



F - 36.596 



E PSX 


EPSY 


K 


A 


1096 


-569 


1 .00 


0. 52 


- 143 1 


V' %<* vU <a» 

O' o' O' 3' 'i' 


-1.31 


4# vC vV 'U 4 # J# 
^ ^ ^ ^ 


-898 


-73 


-0.82 


-11.51 


-616 


-20 


-0.56 


-30. 80 


-25 


420 


-0.02 


0. 06 


-215 


46 8 


-0.20 


0.46 


-514 


525 


-0.47 


0.98 


-312 


861 


-0.28 


0. 36 


-731 


731 


-0.67 


1 .00 


-5 70 


966 


-0.52 


0.59 


-5 70 


966 


-0.52 


0.59 


1 096 


j569 


1 .00 


0. 52 


- 1520 


O' 'fc O' O' O' O' 


-1.39 


^ Jo Jf 

oy' ojs oy* y y op 


-832 


129 


-0.76 


6.45 


- 706 


-106 


-0.64 


-6.66 


-506 


100 


-0.46 


5.06 


-638 


243 


-0.58 


2.63 


-523 


534 


-0.48 


0. 98 


-435 


852 


-0.40 


0. 51 


-483 


988 


-0.44 


0.49 


-629 


1002 


-0.57 


0.63 


-658 


1042 


-0.60 


0.63 



8 

❖ 

0 

0 

3 

5 

0 

6 

5 

1 

1 

8 

•r 

4 

5 

4 

0 

5 

5 

3 

5 

5 
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STRFSS CONCENTRATION DATA 



28 MAY 75 



PLATE # 1205 Y-AXIS 



# OF 


RUNS: 


2 


LOAD 400 


.0 


o 

II 

LL 


. 248 






Y/R 


TNN 


TNO 


NN 


NO 


EPSX 


EPSY 


K 


A 


FF 


-0.5 


-5.0 


51.7 


53.4 


1424 


-676 


l .00 


0.47 


1.0 


-5.8 


-5.0 


63.4 * 


sjcsjc# & # 


-2047 


A J« tl V/ 

-r v 'r v t 


-1.44 


^ -n* nr* -i* 'r* 


1.5 


-0.5 


-5.0 


2.0 


2.1 


227 


328 


0.16 


-0.69 


2.0 


-0.5 


-5.0 


18.5 


22.4 


1 24 


889 


0.09 


-0. 14 


2.5 


-0. 5 


-5.0 


31.3 


28.5 


-537 


74 2 


-0.38 


0.72 


3.0 


-0.5 


-5.0 


35. 2 


34.5 


-489 


947 


-0.34 


0.52 


4. 0 


-0.5 


-5.0 


38.0 


41.5 


-291 


1257 


-0.20 


0.23 


5.0 


-0.5 


-5.0 


45.0 


43.3 


-746 


1034 


-0.52 


0. 69 


6.0 


-0.5 


-5.0 


46.7 


46.2 


-708 


1191 


-0.50 


0.59 


FF 


-0.5 


-5.0 


51.7 


53.4 


1424 


-676 


1.00 


0.47 


1.0 


-0.8 


-5.0 


65.5 * 


%V ^ 

^ *1' 


-1962 


^ X J/ 


-1.38 


a. j. ^ a 

'(• -r ‘i 1 v v 


1.5 


-0.5 


-5.0 


4.6 


4.7 


175 


^380 


0.12 


-0.46 


2.0 


-0.5 


-5.0 


13.0 


22 .4 


1 65 


909 


0.12 


-0.1 8 


2.5 


-0.5 


-5.0 


29. 8 


28.5 


-416 


802 


-0.29 


0.52 


3.0 


-0.5 


-5.0 


35.7 


34.5 


-529 


927 


-0.37 


0. 57 


4.0 


-0.5 


-5.0 


40.3 


41 . 5 


-476 


1165 


-0.33 


0.41 


5.0 


-0.5 


-5.0 


44.2 


43.3 


-682 


1116 


-0.48 


0.61 


6.0 


-0.5 


-5.0 


46.4 


46.2 


-684 


1203 


-0.48 


0.57 



STRESS CONCENTRATION DATA 



29 MAY 7 



PLATE If 1206 

If OF RUNS: 2 LOAD 



Y-AXIS 

200.0 F = 36.596 



Y/R 


TNN 


TNO 


NN 


NO 


E 0 SX 


EPSY 


K 


A 


FF 


-1.0 


-5.5 


55.6 


61.2 


1590 


-481 


1 .00 


0,30 


1.0 


11.2 


-5.5 


55.6 


«JL* U# vC X X 

-y* +1+ -v* 'c 


-1194 


X X X x X X 
V -v V ^ 


-0 . 75 


X •**/ x X 

-v* 


2.0 


-1.0 


-5 . 5 


19.8 


44.7 


1233 


1994 


0.78 


-0.62 


3.0 


-1.0 


-5.5 


3 3.5 


38.0 


- 1 37 


1125 


-0.09 


0. 12 


4. 0 


-1.0 


-5.5 


43.0 


43.1 


-552 


1057 


-0.35 


0.52 


5.0 


- 1.0 


-5.5 


46.0 


49.0 


-443 


1271 


-0.28 


0.35 


6. 0 


-1.0 


-5.5 


48.0 


49.0 


-594 


1198 


-0.37 


0. 50 


FF 


-1.0 


-5 . 5 


55.6 


61.2 


1 5 90 


-481 


1 .00 


0.30 


1.0 


11.2 


-5.5 


55.0 


# 3}: jJ: jfe # 


- 1 200 


X x X X X 

T '1* -V* V* V V 


-0.75 


X V* *** X X 

V ^ -v* 


2.0 


-1.0 


-5.5 


27.7 


34. 7 


106 


1156 


0.07 


-0. 09 


3.0 


-1.0 


-5.5 


36.8 


37.5 


-406 


977 


-0.26 


0.42 


4.0 


-1 .0 


-5.5 


45. 1 


44.5 


- 629 


1057 


-0.40 


0.60 


5.0 


-1.0 


-5.5 


46.6 


44.4 


-744 


997 


- 0 . 4 7 


0.75 


6. 0 


-1.0 


-5.5 


49.0 


47.0 


-777 


1052 


-0.49 


0.74 
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APPENDIX B: COMPUTER PROGRAMS 



B-l MATERIAL PROPERTIES PROGRAM 



EQUIVALENCE (TITLE »RTB(5) 5 

REAL *4 TITLE ( 15) , STRA( 15) , NN , NO , NU, E PS Y < 30) ,SIGMA(30) 
INTEGER*4 P,EPSX 
INTEGERS ITB( 12) /12*0/ 

RFAL*4 RTB( 28)/28*0.0/ 

I TB( 1 ) = I 
I TB( 2) =0 
I TB( 3 ) =6 
I T B ( 4 ) = 6 
I TB( 7 ) = I 
I T B ( 10) =2 
I TB( 1 1 ) = 1 

C INITIAL ALPHANUMERIC DATA 

RE AD (5t 100) STRA 
100 FORMAT! 1 5A4) 

WRITE(6,110) STRA 
110 FORMAT ( *1' » 1 5 A 4 ) 

RE AD ( 5 , 10 0) TITLE 
WRITE (6, 120) TITLE 
120 FORM AT ( 'O' , 15A4) * 

READ! 5,130) K 1 , A , F 
130 FORMAT ( 1 OX , I 2 , 3X , F10 . 5 , 4X , F 1 0. 5 ) 

WRITE (6, 140) K 1 , A , F 

140 FORMAT ( ' 0# OF RUNS:', 12,' AREA »,F6.3,' F = »,F7 

. 3 

DC 30 1=1,10 
READ ( 5 , 100 ) STRA 
30 W R I T E ( 6 , 150) STRA 
READ! 5 , 100) STRA 
WR1TE(6,12 0) STRA. 

RE AD ( 5 , 100) STRA 
WRITE(6, 150) STRA 
150 FORMAT ( » ' ,1 5A4) 

DO 20 J= 1 , K1 



RE AD( 5,155) K2 
155 FORMAT! 5 4 X , I 2 ) 

R EAD( 5,100) STRA 
TARE VALUES 
R EAD( 5 , 160 ) P , TNN , TNO 
160 FORMAT I I 10,2F10.4) 

I S IGM A= P/ A 

WRITE(6,170) P , I SIGMA , TNN, TNO 
170 FORMAT l '0',14,1X,I5,2(1X,F6.1)) 

DATA REDUCTION 

DO 10 I = 1 , K2 

READ ( 5 , 160 ) P , NN , NO 
I SIGMA=P/A 

S I GM A ( I ) = F LO AT { ISIGMA) 

DELNN=NN-TNN 
DEl.NC= NO- T NO 

1 EPSY=F*( 1 o5*DELNQ-DELNN) 

E P SY ( I ) =FLQAT( IE PSY) 

E PSX = F * ( 1 . 5*DELN0- 2 . *DELNN ) 

NU=-EPSX/EPSY( I ) 

10 W R I T E ( 6 , 1 8 0 ) P, I SIGMA, NN, NO, DELNN,DELNO, I E P SY , E t> 5 X , NU 
1 BO F0P3‘AT(' ' , I 4, 1 X , I 5,4( IX, = 6. 1 ) , 2 ( 2X, I ) , 2X, F4.3) 

IF( J.EQ.K1 ) WRITE <6 i 110) 

C PLOT STRESS-STRAIN CJRVE 

IF(J.GT.l) I T 3 ( 1 ) = 2 
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IF(J.EQ.Kl) I T8 ( 1 ) =3 
I TB ( 2 ) = I TB ( 2 ) + 1 
IF < ITB< 2 ) . E0.5 ) I TB ( 2 ) = 1 
CALL D RAW P(K2,EPSY, SIGMA, I T B , RT B ) 
20 CONTINUE 
STOP 
END 



) 



59 



1 00 

no 

120 

130 

140 

30 

150 

155 

160 

1 70 

165 

6 



B-2 STRESS CONCENTRATION PROGRAM 



EQUIVALENCE ( T I T L E , RT B ( 5 ) ) 

REAL*4 TI TLE ( 15) , STR A( 1 5) , NN,NO, K ( 30 ) , X < 30 ) , A( 30) ,KF 
INTEGER *4 EPSX,FEPSY 
INTEGERS ITBC 121/12*0/ 

RE AL*4 RTB(28)/28*0.0/ 

ITBC 1 ) = 1 
I TB ( 2 ) = 3 
I TB ( 3 ) = 6 
I TB( 4) =6 
I TB( 7 ) = 1 
I TB( 10) =2 
ITBC 11 ) = 1 
DO 1 ICCUNT=1,3 

INITIAL ALPHANUMERIC DATA 
READ! 5 » 100) STRA 
FORMAT! 15A4) 

WRITE(6»110) STRA 
FORMAT! « 1 • , 15A4) 

R E AD ( 5 » 1 0 0 ) TITLE 
WRITE (6» 120) TITLE 
FORMAT! 1 G ' t 1 5A4) 

READ! 5,130) K1 , P ,F 

FORMAT! 10X,I2,SX,F10.5,4X,F10.5) 

WRITE! 6, 140) K1,P,F 

FORMAT! '0# OF RUNS:', 12,' LOAD ',F6.1,’ F= ',F7. 

3) 

DC 30 1=1,10 
READ! 5, 100) STRA 
WRITE (6, 150) STRA 
READ! 5 , 100 ) STRA 
WRITE!6,120) STRA 
RE AD (5, ICO) STRA 
WRITE (6, 150) STRA 
FORMAT!' *,15A4) 

DO 20 J=1,K1 



READ! 5, 155) K2 
FORMAT! 54X , 12) 

READ! 5 , 100 ) STRA 

FAR FIELD VALUES 
RE AD (5, 160) TNN , TNO , NN , NO 
FORMAT ( 10X,4F10. 4) 

DELNN= NN-TNN 
DELNO= NO- TNO 

FEPSY=F*( 1 ,5*DEL NO-DEL NN ) 

E P S X = F * ( 1 .5*DELN0-2.*DELNN) 

KF=1 . 0 

AF=-EPSX/ FLOAT ( FEPSY ) 

WRITE! 6 , 1 70) TNN , TNO, NN, NO, FEPSY , E PSX , <F , AF 
FORMAT! 'OFF ',2(F5.1,1X),2(F6.1,1X),2!I6,1X),F5.2,1X, 
F6 . 2 ) 

DATA REDUCTION 
DO 10 1=1, K2 

RE AD { 5 , 16 5) X ( I ) , TNN , T NO , NN , NO 
FORMAT! 5 FI 0.4) 

DELNN= NN-TNN 
DE LNO= NO- TNO 

I E PSY=F* ( 1 . 5*DEL NO-DELNN ) 

E PSX= F* ( 1 . 5*DELN0-2.*DELNN) 

K( I ) =IEPSY/FLOAT (FEPSY) 

A ( I ) =-EPSX/ FLOAT ( I FPSY ) 

IF(NO.NE.O.O) GO TO ID 
I E PSY = F *D E LMN/ 1 .36 
N0= 10**8 
EPS X =10* *8 

K ( I)=IEPSY /FLOAT (FEPSY) 
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A ( l)=10**8 

10 WRI TE < 6 i 1 30 ) X( I ) , INN, TNO, NN,NO, I EPSY , E PSX , K( I ) , A( I ) 
180 FORMAT ( » • ,F3. 1 , IX ,2( F8. 1 , IX) , 2 ( F6. 1 » IX) , 2 ( 16, IX) ,F5.2 

t IX t F6.2 ) 

IF(J.EQ.Kl) WRITE(6,110) 

C PLOT STRESS-STRAIN CJPVE 

IFIJ.GT .1 ) I TS ( 1 ) = 2 
IF ( ICOUNT.NE .3 ) GO TO 2 
IF(J.EQ.Kl) I TB ( 1 ) =3 
2 CONTINUE 

CALL DRAWP(K2,X,K,ITB,RTB) 

20 CONTINUE 

ITB(2)=ITB(2)+1 
I F < I T B < 2) . EQ . 6 ) I TB C 2 ) = 1 
1 CONTINUE 
STOP 
END 
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